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We extract the Q-Q potential by using the thermodynamic quantities obtained in lattice gauge 
calculations. The potential is tested and found to give spontaneous dissociation temperatures that 
agree well with those from lattice gauge spectral function analysis. Using such a Q-Q potential, we 
examine the quarkonium states in a quark-gluon plasma and determine the 'quark drip lines' which 
separate the region of bound color-singlet QQ states from the unbound region. The characteristics 
of the quark drip lines severely limit the region of possible bound QQ states with light quarks to 
temperatures close to the phase transition temperature. Bound quarkonia with light quarks may 
exist very near the phase transition temperature if their effective quark mass is of the order of 
300-400 MeV and higher. 
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I. INTRODUCTION 



. The degree to which the constituents of a quark-gluon plasma (QGP) can combine to form composite entities 

■ is an important property of the plasma. It has significant implications on the nature of the phase transition, the 
, quark-gluon plasma equation of state, the probability of recombination of plasma constituents prior to the phase 

PsJ ■ transition, and the chemical yields of the observed bound hadrons. The successes of the recombination model P, 
i 0, H, 0, H, H, 0, H, B suggests that quark partons may form bound or quasi-bound states in the quark-gluon 

■ plasma, at least at, or close to, the phase transition temperature. It is an important theoretical question as to the 
' range of temperatures in which thes e q uarkonia may be bound or quasi-bound. The successes of the thermal model 

^ ■ Q, 1, 1 i M M M M M, M M,MMM for hadron yields also raises the important question whether hadrons 
' may become bound or quasi-bound in the quark-gluon plasma. If they arc indeed bound in the quark-gluon plasma, 
Oh, the approach to chemical equilibrium may commence in the quark-gluon plasma phase before the phase transition 
O • and the boundaries of the quark-gluon plasma phase and the hadron phase may overlap. 

, Recent spectral analyses of quarkonium correlators indicated that J/ip may be stable up to l.GTc where Tc is 
• the phase transition temperature (2ll . [2^ . Subsequently, there has been renewed interest in quarkonium states in 
i> ! quark-gluon plasma as Zahed and Shuryak suggested that QQ states with light quarks may be bound up to a few Tc 
[23|. Quarkonium bound states and instanton molecules in the quark-gluon plasma have been considered by Brown, 
Lee, Rho, and Shuryak [24| . As heavy quarkonia may be used as a diagnostic to ol (251, there have been many recent 
investigations on the stability of heavy quarkonia in the plasma [2^, [26 .[27l [28l . [29l . I30l.l3ll. [3^ . [ssl . [U, [H, H^, [13, [H, [1^ . 

Previously, DeTar [i^l, Hansson, Lee, and Zahed [4l|, and Simonov 42l| observed that the range of strong interaction 
is not likely to change drastically across the phase transition and suggested the possible existence of relatively narrow 
low-lying QQ states in the plasma. On the other hand, Hatsuda and Kunihiro [1^ considered the persistence of 
soft modes in the plasma which may manifest themselves as pion-likc and sigma-like states. The use of baryon- 
strangeness correlation and charge fluctuation to study the abundance of light quarkonium states in the plasma have 
been suggested recently [i^ . [45j . 

We would like to investigate the composite properties of the plasma and to determine its 'quark drip lines'. We shall 
focus our attention on color-singlet Q-Q states and the Q-Q quantities in this paper refer to those of color-singlet Q-Q 
states unless specified otherwise. Here we follow Werner and Wheeler [4^ and use the term 'drip line' to separate 
the region of bound color-singlet QQ states from the unbound region of spontaneous quarkonium dissociation. It 
should be emphasized that a quarkonium can be dissociated by collision with constituent particles to lead to the 
corresponding 'particle-dissociation lines', which can be interesting subjects for future investigation. 

We would like to use the potential model to study the stability of quarkonia, as the potential model can be used 
to evaluate many more quantities than the lattice gauge spectral function analysis. The potential model lends itself 
to extrapolation into unknown regions of quark masses and temperatures. An important physical quantity in the 
potential model is the Q-Q potential between the quark Q and the antiquark Q at a se par ation i? at a temperature T. 
Previous work in the potential model uses the color-singlet free energy Fi{R, T) [26l. [Sll. |38| or the color-singlet internal 
energy Ui{R,T) [22,,^, i37j obtained in lattice gauge calculations as the Q-Q potential, without rigorous theoretical 
justifications. Here, the subscripts of Ui{R,T) and Fi{R,T) refer to the color-singlet property of the Q-Q system. 
The internal energy Ui{R,T) is significantly deeper and spatially more extended than the free energy Fi{R,T). The 
degree of quarkonium binding will be significantly different whether one uses the internal energy Ui{R,T) or the 
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free energy Fi {R, T) as the Q-Q potential. Treating the internal energy Ui {R, T) as the Q-Q potential led Shuryak 
and Zahed to suggest the possibility of color-singlet quarkonium states with light quarks in the plasma [23| . The 
conclusions will be quite different if one uses the free energy Fi{R,T) as the Q-Q potential. 

While Fi{R,T) or Ui{R,T) can both be used as the Q-Q potential at T = (at which Fi{R,T) = Ui{R,T)), the 
situation is not so clear in a thermalized quark-gluon plasma. It is important to find out the meaning of these ther- 
modynamical quantities calculated in the finite-temperature lattice gauge theory so as to extract the Q-Q potential. 

If one constructs the Schrodinger equation for the Q and Q in a thermal medium, the Q-Q potential in the 
Schrodinger equation contains those interactions that act on Q and Q, when the medium particles have re-arranged 
themselves self-consistently. On the other hand, the total internal energy Ui{R, T) contains not only these interactions 
that act on Q and Q, but also the internal gluon energy Ug{R,T) relative to the gluon internal energy UgQ in the 
absence of Q and Q, as shown deductively in Ref. [l^l starting from the definition of the free energy in lattice gauge 
theory in quenched QCD. If the gluon internal energy Ug{R,T) were independent of R, then Ui{R) could well be 
used as the Q-Q potential. However, in the grand canonical ensemble, Ug{R,T) depends on R. To get the Q-Q 
potential, it is therefore necessary to subtract Ug{R,T) — Ugo from Ui{R,T). As the subtleties of these results may 
not appear evident and the problem of non-perturbative QCD so intrinsically complicated, a thorough understanding 
of an analogous, but not identical, problem in QED is worth having. Therefore, we examine in detail the simple QED 
case of Debye screening of charges Q and Q in a massless charged medium in a grand canonical ensemble, where 
the results can be readily obtained analytically. We would like to show that there is a relationship between the Q-Q 
potential and the total internal energy when screening occurs; the Debye screening Q-Q potential between two static 
opposite charges in QED can be obtained from the total internal energy by subtracting out the internal energy of the 
medium particles. 

The results in the Debye screening case in QED support our previous conclusion in Ref. [l^ that in the QCD 
lattice gauge calculations in the grand canonical ensemble, it is necessary to subtract out the i?-dependent internal 
energy of the QGP from the total internal energy in order to obtain the potential between Q and Q in the plasma. 
Additional lattice gauge calculations may be needed to evaluate the QGP internal energy in the presence of Q and Q. 
It is nonetheless useful at this stage to suggest approximate ways to evaluate the QGP internal energy. We proposed 
earlier a method by making use of the equation of state of the quark-gluon plasma obtained in an independent lattice 
gauge calculation [231. The equation of state provides a relationship between the QGP internal energy and the QGP 
entropy content. As the QGP entropy content is the difference Ui — Fi, the Q-Q potential can be represented as 
a linear combination of Ui and Fi, with coefficients depending on the quark-gluon plasma equation of state. The 
proposed potential was tested and found to give spontaneous dissociation temperatures that agree well with those 
from lattice gauge spectral function analysis in the quenched approximation. 

The comparison of the potential model results with those from spectral analyses in the same quenched approximation 
is useful as a theoretical test of the potential model. However, in the quenched approximation, the quark-gluon plasma 
is assumed to consist of gluons only and the effects of dynamical quarks have not been included. As dynamical quarks 
provide additional screening, one wishes to know whether this additional screening will modify the binding energies of 
quarkonia significantly or not. The presence of dynamical quarks also lowers the phase transition temperature from 
269 MeV for quenched QCD to 154 MeV for full QCD with three flavors [l^. For these reasons, it is necessary to 
include dynamical quarks to assess their effects on the stability of quarkonia. The knowledge of the single-particle 
states using potentials extracted from lattice gauge calculations in full QCD can then be used to examine the stability 
of both heavy and light quarkonia and to determine the location of the quark drip lines. We focus our attention 
mainly on heavy quarkonia for which a non-relativistic treatment is a good description. However, the problem of 
the stability of quarkonia with light quarks is intrinsically so complicated and the question of their stability up to 
a temperature of few units of Tc so important [H-fl^ that even an approximate estimate using the non-relativistic 
potential model is worth having. The subject of light quarkonia will be examined again, with the inclusion of the 
relativistic effects as in recent works (48l.l49j. in the course of time. 

The authors of Refs. [13, [5l| claim recently that potential models cannot describe heavy quarkonia above T^, as 
their potential model correlators fail to reproduce lattice gauge correlators for all cases with all types of potentials. 
Such a complete disagreement for all cases and all types of potentials suggests that the lack of agreement may not 
be due to the potential model (or models) themselves but to their method of evaluating the meson correlators in 
the potential model. We show recently that when the contributions from the bound states and continuum states are 
properly treated [lOl , the potential model correlators obtained with the proposed potential in Ref. [23| are consistent 
with lattice gauge correlators. 

This paper is organized as follows. In Section II, we describe the puzzling behavior of the increase of the entropy of 
the QCD medium as the separation between Q and Q increases. To understand such a behavior, we introduce in Section 
HI a simple model of Debye screening in QED for which various thermodynamic quantities can be readily calculated. 
In Section VI, the variation of the number density and entropy density of the medium particles for a Q-Q pair in Debye 
screening is shown to depend on the Q-Q separation R when the second-order contributions are included. In Section 
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V, we find similarly that the internal energy of the medium particles in Debye screening also increases with R when 
second-order contributions are included. In Section VI, we examine the Schrodinger equation for the relative motion 
of Q and Q and identify the potential between Q and Q in Debye screening. We reach the conclusion that in order 
to obtain the Debye screening potential between two static charges in the grand canonical ensemble, it is necessary 
to subtract out the internal energy of the medium particles from the total internal energy. Returning to lattice 
gauge calculations in Section VII, we suggest a method in which the medium internal energy can be approximately 
determined and subtracted by using the information on the quark-giuon plasma of state. Consequently, the Q-Q 
potential turns out to be a linear combination of Ui{R,T) and T), with coefficients depending on the equation 

of state. In Section VIII, we use different potentials to calculate the spontaneous dissociation temperatures for various 
quarkonia and compare them with those obtained from lattice gauge spectral function analyses in quenched QCD. 
The dissociation temperatures obtained with the proposed linear combination of C/i and Fi give the best agreement 
with those from the spectral function analyses. In Section IX, we show how the thermodynamical quantities in full 
QCD with two flavors are parametrized. The dissociation temperatures for various heavy quarkonia for full QCD 
with two flavors are obtained in Section X. We introduce the quark drip lines in quark-giuon plasma in Section XI. 
We conclude and summarize our discussions in Section XII. 



II. THERMODYNAMICAL QUANTITIES IN LATTICE GAUGE CALCULATIONS AND DEBYE 

SCREENING 

Thermodynamical quantities for a heavy quark pair in the color-singiet state was studied by Kaczmarek and Zantow 
in quenched QCD and in full QCD with two flavors [H, [s^l- They calculated (tri(r/2)it(— R/2)) and obtained the 
color-singlet free energy i^i(R, T) from 

(irL(R/2)L+(-R/2)) = e-^i(i^>^)/'=^. (1) 

Here trL(R/2)L^(— R/2) is the trace of the product of two Polyakov lines at R/2 and — R/2. The free energy 
i^i(R, T), in the presence of the Q-Q pair, is measured relative to the free energy without the Q-Q pair. The quark 
and the antiquark lines do not, in general, form a close loop. As a gauge transformation introduces phase factors 
at the beginning and the end of an open Polyakov line, (trL(R/2)it(— R/2)) is not gauge invariant under a gauge 
transformation. Calculations have been carried out in the Coulomb gauge which is the proper gauge to study bound 
states. 

From the free energy Fi, Kaczmarek and Zantow [H, [s^l calculated the internal energy Ui using the statistical 
identity 

UiiR,T)=Fi{R,T)+TSi{R,T), (2) 

where Si{R,T) = —dFi{R,T)/dT is the entropy of the system in the presence of a color-singlet Q-Q pair and is 
measured relative to the entropy of the system in the absence of the Q-Q pair. 

In order to extract the Q-Q potential from thermodynamical quantities calculated in the lattice gauge theory, we 
need to understand the behavior of the free energy Fi{R,T), the internal energy Ui{R,T), and T times the entropy, 
TSi{R,T), which we shall also abbreviatingly call "the entropy" for simplicity of nomenclature. As these three 
quantities are related by Eq. ^ we need to find out the behavior of only two of these three quantities. Following the 



terminology used in lattice gauge calculations |34| . Ui, Fi, and TSi are defined as being measured relative to their 
corresponding quantities in the absence of Q and Q. 

We can begin by studying the entropy TSi {R, T) of the system. We note that the lattice gauge calculations show 
that the total entropy TSi {R) increases as a function of i?, and saturates after the separation reaches a large value 



of R [3JI as shown in Fig. 1(a). What does such a behavior tell us about the response of the medium particles to the 
presence of the external color sources of Q and Q in a thermal bath? 

In the system under consideration, the system consists of the Q, the antiquark Q and the quark-giuon plasma. For 
simplicity we can examine the quenched case for which the quark-giuon plasma is assumed to consist of giuons only. 
The entropy of the system therefore comes from the sum of the entropies of Q, Q and the giuons. However, in the 
lattice gauge calculations in a thermal bath, the Q and the Q are held fixed and do not contribute to the entropy of 
the system. The entropy of the system TSi comes entirely from the giuons. By fixing a temperature and focusing 
our attention at the state of thermal equilibrium in lattice gauge calculations, the giuons are in a grand canonical 
ensemble in contact with the thermal bath. The content and thermodynamical properties of the giuons are determined 
by the condition of thermal equilibrium at the fixed temperature of the thermal bath. The observed behavior of the 
entropy TSi{R) as a function of R in Fig. 1(a) suggests that the gluon entropy content increases as the separation R 
increases, and the entropy saturates when R reaches a certain limit. 
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FIG. 1: (a) The total entropy TSi{R,T) as a function of i? at T = 1.3Tc from the lattice gauge calculations [H. (6) The 
ratio ri{R) = TS{R)/TS{R oo) (which is also equal to N{R)/N{R ^ oo) and U{R)/U{R oo)) as a function of ii/r^^ in a 
thermal medium under Debye screening. 



As wc know from the work of Landau and Belcnkii 52] , the entropy content of a system is closely correlated with 
the number content of the particles in the system. The behavior of the entropy suggests that the number of gluons 
increases as the separation R increases. How do we understand such a behavior? If the number of gluons increases 
as a function of R, what happens to the internal energy content of the gluons as a function of R, and how does that 
affect the total internal energy and the potential between Q and Q we wish to extract? 



III. ANALOGOUS PROBLEM OF DEBYE SCREENING IN QED 



The QCD problem of a quark and an antiquark in the presence of the actions of the gluons is so complicated that 
it is worth having even a good understanding of the puzzling behavior in an analogous, but not identical, problem 
of a static positive charge Q and another negative charge Q under the action of charged medium particles in QED. 
We would like to ask whether there is a similar behavior of the number and the entropy of medium particles as the 
separation R between the positive charge Q and the negative charge Q increases. If the number of medium particles 
increases as a function of R, what happens to the internal energy content of the system and the medium particles as 
a function of R, and what is the relationship between the potential between Q and Q and total internal energy we 
wish to extract? 

Accordingly, we study the simple system of Debye screening and consider a Q with charge +q held fixed at — R/2 
and a Q with a charge —q at R/2, in a medium of massless fermions with charges e± = ±q, in a thermal bath of 
temperature T in a grand canonical ensemble. Any pair of charged particles with charges ei and 62 separated by a 
distance r are assumed to interact with a Coulomb interaction eie2/|r|. To make the problem simple, we shall assume 
the attainment of local thermal equilibrium. 

When the external charges Q and Q are introduced into the medium, the medium fermions will re-arrange themselves 
in both coordinates and momenta to reach a new local thermal equilibrium. The total number of the medium particles 
in the system is not a constant of motion but is determined by the condition of thermal equilibrium, maintained by 
the thermal bath. 

When the external charges Q and Q are separated by a distance R at a temperature T, the self-consistent re- 
arrangement of the medium charged fermion particles leads to a potential I^(r,R) at a point r and under a local 
thermal equilibrium, the momentum distribution of the medium particles at r in the Born-Oppenheimer treatment is 
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given by 



/±(r,p,R) = j^— (3) 

exp{[p + e±y(r,R)]/T} + 1 



Here and henceforth, the + and — subscripts designate quantities for the positive and negative medium particles 
respectively. 

From the above Wigner function distribution, we can obtain various thermodynamic quantities. The integration 
of the Wigner function over all momenta gives the spatial number density distribution n±(r, R) at r, when Q and Q 
are separated by R, 

n±{r,R) ^ ^ J dp f±{r,p,IL), (4) 

where g is the degeneracy of the levels. We can consider the high temperature case for which it is useful to expand 
various quantities as a power series of e±V{{r,'R.)/T. Up to the second order in [e±y(r, R)/T]^, the medium particle 
density is 

n±(r,R) =nl{l- ai [e±V^(r, R)/r] + a2[e±y(r, R)/r]2} , (5) 

where 

4-|^4C(3)r(3), (6) 



K(2) 
|C(3) 



ai = 1^ = 0.91233, (7) 



and 



ic(i) 

4C(3) 



a2 = 1^ = 0.3845. (8) 



In passing, we note that if the medium particles obeys Boltzmann statistics, the coefficients oi and 02 would equal 
oi = 1 and a2 = 0.5, as it follows from the expansion of the well-known Boltzmann distribution n-|-(r,R) = 
rij. exp{— e±y(r, R)/T} for Boltzmann particles in an external field. The values of the ai and 02 coefficients for 
Boltzmann statistics differ only slightly from the corresponding values in ([7]) and ([5]) for fermion particles. 

For the Fermi-Dirac medium particles, the entropy density at r, when Q and Q are separated by R, is given by 

a±(r,R) = 1^ jp" dp{-f± ln/± - (1 - /±) ln(l - /±)} . (9) 

Upon substituting the Fermi-Dirac distribution of Eq. ([3]) into the above equation, we find 

a±(r,R) = ^ J p' dpf±{r,p,R) |^ + e±y(r,R)| ^ T 

^ (|+e±y(r,R))/r, (10) 



and we obtain 



where 



a±(r,R) =al{l- 6i[e±y((r, R)/T] + b2[e±V{r,R)/T]^} , (11) 



,0 _ .9^3 7,,, ,4, 



al = ^-C(4)-r(4), (12) 
± 27r2 8^^ ^3 ^ ^' ^ ^ 



|C(3)[|r(4)-r(3)] , , 

61 = y ' V ;j ^ 0.7139, 13 



C(4)|r(4) 
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and 

'^^ ^ ^ = 0.2171. (14) 



K(4)|r(4) 

Here we have purposely written the numerators of the h\ and 62 coefficients as a difference where the first term comes 
from {'^p/'i) and the second term comes from {e±V) of Eq. (fTO|l . 

The evaluation of various thermodynamical quantities requires the knowledge of T^(r,R). To determine y(r,R) 
self-consistently, we have the charge density at the point r, when Q and Q are separated by R, 

R R 

ptotai(r, R) = q6{Y + y) - q5{r - y) + e+n+(r, R) + e_7i_ (r, R). (15) 

Using the number density distribution given in Eq. ([5]), the charge density, up to the second power in e±y(r,R)/T, 
becomes 

R R 

Ptotai(r,R) = q6{v + qS{r - y) - noaiq'^V{r,R), (16) 

where uq = + nP_ and the zeroth-order and second-order terms of Eq. (O cancel out on account of the presumed 
charge neutrality of the system for which = rfi_. The Poisson equation for the potential is then given by 

VlV{v, R) = -4^ |g 5{v - q 5{v - n^a.q^Viv, R)| , (17) 

which has the solution 

y(r,R) = ^ ^ , (18a) 

r_|_ r_ 



r± = |r±R/2|, (18b) 
/47rnoaig2 l 

/^ = V^— = 

where fi is the Dcbye mass and is the Dcbye screening length. 

IV. VARIATION OF NUMBER DENSITY AND ENTROPY DENSITY WITH R IN DEBYE 

SCREENING 

The simple solution F(r,R) in Eq. ([T5)l of the last section allows us to have a profile of the self-consistent medium 
particle number density and entropy density in all spatial points at local thermal equilibrium at T. In Eqs. ^ and 
(|lip . the coefficients of oi, 02, 61 and 62 are all positive. For positive medium particles, the first-order increment in the 
number density n+(r,R) and CT-|-(r,R), [also M+(r,R) in Eq. ((32)) ] are therefore measured by [— e+V^(r, R)/T] illus- 
trated in Fig. 2 as a function of p/r^ and z/r^^ . As one observes, the first-order contributions, given by — e+y (r, R)/T, 
represent a depletion for the positive medium particles near the positive static charge Q at — R/2, and an enhance- 
ment near the negative static charge Q at R/2. The degree of depletion and the degree of enhancement are equal and 
opposite to each other. When we sum over all spatial points, the sum of the first-order depletion and enhancement 
cancel each other to give a zero total contribution. 

The second-order contributions are proportional to [e+y(r, R)/T]^ and are always positive. They are illustrated 
in Fig. 3 as a function of p/r^^ and z/r^^ for different Q-Q separations. They always enhance the number density and 
the entropy density. The enhancement is small when the two static charges are close together in Fig. 3(a), as there is 
a substantial cancellation of the two terms in Eq. p8b ). The enhancement reaches a constant value when the static 
charges Q and Q reaches a separation of 1-2 units of the Dcbye screening length as shown in Fig. 3b and 3c. 

If one integrates over all spatial points to obtain the total number of positive charge medium particles, one finds that 
the integration over the first-order term, J dr[— e+V^(r, R)/T], is zero because the depletion cancels the enhancement. 
However, the second-order contributions always give a positive contribution, and the number of positive medium 
particles, measured relative to its corresponding quantity in the absence of Q and Q, is given by 

N+{R)^ J dr{n+{r,TV) -nl) =nl J dr a2[e+V{r,R)/Tf . (19) 
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(a) R = 0.1r^, (b) R = l.Or^, (c) R=1.9r^ 

FIG. 2: The first-order term [—e+V{pz,T)/T] in units of q'^/r^ in Eqs. (O, and p2|l that contributes to the increments 

in n+{pz, R),a+{pz, R), and u+{pz, R), as a function of tire spatial coordinates p/r^^ and a/r^. Fig. 2(a) is for R = O.lr^, Fig. 
2(6) for R = l.Or^ , and Fig. 2(c) for R = 1.9r„ . 



We can write the above as 

N+{R) 



N+{R ^oo) 

where 



r,iR) = l-JdC 



v{R) (20) 



(21) 



C+ C- 



with C = r/^^D, C± |i"± R-/2|/rD, ?7(0) = 0, and ri(R ^ oo) = 1. In Eq. ([20]), ^+(i? ^ oo) is the increment in the 
number of positive medium particles when Q and Q are far separated, 



N+[R^oo)=nl3a,(^^^ (22) 

The increase in the total number of positive medium particles is small when Q and Q are close together, and the 
increase saturates when R reaches a few units of the Debye screening length . For our charge-neutral system, we 
obtained from Eq. Q in a similar way N^{R) = N-^{R). 

Similarly, the entropy density of the medium particles is depleted near the static charge of the same sign, and is 
enhanced in the vicinity of the static charge of the opposite sign. When integrated over all spatial points, the depletion 
and the enhancement cancel to the first order. The second-order contributions always give a positive total entropy, 

TS±iR)=T J dr{a±(r,R)-4} = J dr b2[e±Vir,R)/T]^ , (23) 

which is small at small R and saturates at large R. We can write the above as 

TS±{R) 



TS±{R^ oo) 

where 



= viR) (24) 



TS^{R^^)=Tal3b,(-^] ^r^ (25) 
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(a) i? = 0.1r^, (b) i?=1.0r^, (c) R=1.9r^ 

FIG. 3: The second-order term {e+V{pz, R)/Tf in units of {q'^/r^Tf in Eqs. (O, (HU, and ((32]) that contributes to n+{pz,R), 
a+{pz, R), and u+{pz, R), as a function of the spatial coordinates p/r^ and z/r^-,. Fig. 3(a) is for R — O-lr^^, Fig. 3(fe) for 
R = l.Or^ , and Fig. 3(c) for R = 1.9r^ . 



A comparison of N± and TS± shows that 

N±{R) TS±{R) 



N±{R ^oo) TS±{R -^oo) 



viR)- (26) 



If wc define iV(i?) = N+ + N-{R), and TS{R) = TS+{R) + TS-{R), the increment of the total number and entropy 
of the medium particles due to the presence of Q and Q (measured relative to the corresponding quantities in the 
absence of Q and Q), is 

- -,iR). (27) 



N{R^ oo) TSiR->oo) 

Thus, the ratios TS{R)/TS{R oo), N±{R)/N±{R oo), and TS{R)/TS{R oo) [also U{R)/U{R oo) as we 
shall see in the next section] behave in the same way as a function of R. We show the behavior of ijiR) = 
TS{R)/TS{R oo) for the Debye screening case in Fig 1(6) and it has the same shape as TSi{R) obtained in 
the lattice gauge calculations shown in Fig. 1(a). 

One can therefore understand that as the result of constraining the system to be in contact with a thermal bath 
in the grand canonical ensemble, the medium particle numbers and entropy increase with increasing R to maintain 
a thermal equilibrium until they saturate at large separation. The thermal bath is therefore a participant in altering 
the content of the medium particles, when the Q is separating from the antiquark Q. (See Section VI for a discussion 
on the role of the thermal bath). 



V. VARIATION OF THE INTERNAL ENERGY WITH R IN DEBYE SCREENING 

Previously in quenched QCD, we find a relationship between the total internal energy and the Q-Q potential [131 . As 
the problem of non-perturbative QCD so intrinsically complicated, a thorough understanding of such a relationship 
in an analogous, but not identical, problem in QED is worth having. We are therefore motivated to examine the 
mechanism of Debye screening due to the medium particles in QED. Of particular interest is to see whether the 
relationship between the total internal energy and the Q-Q potential in Debye screening in QED resembles a similar 
relationship in lattice gauge theory obtained previously in QCD [27l |. 



9 



When Q and Q separated by a distance R are screened by medium particles, the total internal energy t/totai(^) 
of the system of medium particles, Q, and Q is the sum of the kinetic energy of the medium particles and the 
interaction energies of the the medium particles, Q, and Q, when the medium particles have re- arranged themselves 
self-consistently, 

UtoUR) = j '^""{■ir^j P'dp[.f+ir, P, R) + /- (r, P, R)]p| 

+ i y" dr|e+<5(r+5:) + e_5(r-^) + e+n+(r,R) + e_n_(r,R)|F(r,R) (28) 
dr{ul + u"_), 



where u^. is the internal energy density of the medium in the absence of Q and Q given by 

4 = |5^C(4)r(4), (29) 

and ?7totai(^) is measured relative to the total internal energy of the system in the absence of Q and Q, C^^odium = 
/ (ir(u^j_ + u^). Using the solution of V{r, R) in Eq. ((TS]) and excluding the infinite energies of a point source acting 
on itself, Eq. (^5]) gives 



t/totai(i?) = I '^'^ [{^+(r,R) - <} + {^-(r, R) - ^° }] , 



(30) 



where u±(r,R) are the internal energy density of the positive and negative charged medium particles in the presence 
of Q and Q given by 



u±(r,R) = ^ J dpf±{r,p,R) l^p 



e±V{v,Il) 



= {P+ 2 ^ ' 

Using the local Fcrmi-Dirac distribution of Eq. ([3]) , we obtain 

u±(r,R) = ul {1 ~ ci[e±V^((r,R)/r] +C2[e±F(r,R)/r]2} , (32) 
where ci and C2 are positive constants, 

ja3)|r(4)-ffili 
iC(4)r(4) 

and 

^2 = |C(4) ' = ^-^^^^^ 

A comparison of u±(r,R) in Eq. ([32|) with n±(r, R) and u±(r,R) in Eqs. ^ and (fTTj) indicates that the first-order 
and second-order contributions to u±(r, R) behave in the same way as those of n±{r, R) and (T±(r, R). The first-order 
increment of the internal energy density of the medium particles is suppressed near the static charge of the same sign, 
and is enhanced near the static charge of the opposite sign as shown in Fig. 2, while the second-order contributions 
to the internal energy are always positive as shown in Fig. 3. The first-order contributions cancel each other, when 
they are integrated over all the spatial points. However, the second-order contributions are always positive, and the 
integration of the second-order contributions always yield a positive quantity. The total medium internal energy of 
the positive and negative medium particles, measured relative to the corresponding quantities in the absence of Q 
and Q, are 

U±iR)^ J dr{u±{r,'R)-ul} =nl J dr C2[e±F(r, R)/r]2. (35) 
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Relative to the total internal energy in the absence of Q and Q, the total internal energy in the presence of Q and Q, 
is 

f/medium(i?) " f/°cdium = U+iR) + C/_(i?). (36) 

Eq. pO|) can therefore be written as 

Utotal{R) = 1- C/mcdmm(^) — C^modium (37) 

From Eqs. ((35)) . p9|) . and pT|) . the quantity U±{R) can be written as 

C/±(i?) 



C/±(i?^ cx)) 

where 



viR), (38) 



[/±(i?^c^)=43c2(^') yr^. (39) 



A comparison of N± and TS± shows that 

U±{R) _ N±{R) _ TS±{R) 



U±iR ^oo) N±{R ^oo) TS±{R ^oo) 

Consequently, if we define U{R) = U+{R) + U-{R), we also have 

U{R) _ N{R) _ TSjR) _ 
U{R ^ oo) ^ N{R ^ oo) ^ TS{R ^ oo) ~ 



7j{R). (40) 



V{R)- (41) 



In this simple model of Debye screening, the ratios of TS{R)/TS{R oo), N{R)/N{R -> oo) and U{R)/U{R oo) 
are equal and their behavior is shown in Fig. 1(6). The entropy, total number, and internal energy of the medium 
particles (relative to the corresponding quantities in the absence of Q and Q) is zero at i? = and increases as a 
function of R until they saturate when the separation R reaches a few units of the Debye screening length. 



VI. SCHRODINGER EQUATION FOR Q-Q SYSTEM IN DEBYE SCREENING 

As the dynamics of a quark and an antiquark in a quark-gluon plasma is very complicated, it is worth having a good 
understanding of the mechanism of screening and its effects on the interaction between a Q and an antiquark Q in an 
analogous, but not identical, problem. We are therefore motivated to examine the mechanism of Debye screening due 
to the medium particles in QED. Of particular interest is to obtain the relationship between the Q-Q potential in the 
Hamiltonian and the total internal energy J/totai, in order to find out whether this relationship resembles a similar 
relationship in lattice gauge theory obtained previously in Ref . [l^l . 

It is however somewhat tricky to determine the Hamiltonian for the Q and Q system in a medium under Debye 
screening (and analogously, but not identically, under color-charge screening in QCD). We can follow the basic 
principles of statistical physics as described by Landau and Lifshitz [ssj . Accordingly, we start with a closed system 
of a heavy Q and Q with medium particles and consider a small "subsystem" S that contains the Q, the Q, and the 
medium. As the number of medium particles of the whole closed system is very large, the number of medium particles 
contained in the small subsystem S can still be very large, and a statistical description of this small subsystem S is 
applicable. This subsystem S is not a closed system and it undergoes all kinds of interaction and medium particle 
exchanges with the complementary part S' of the whole system. We can describe this subsystem S to be in contact 
with a very large complementary part 5", which we can call a "thermal bath" in this connection. The QED system 
discussed in the last few sections or the QCD system in lattice gauge calculations, (the so-called "system" of Q, Q, 
and medium particles in contact with a thermal bath) , corresponds in actual fact to the "subsystem" S out of the 
whole closed system S + S' . 

Under a perturbation of the subsystem S away from thermal equilibrium such as occurs in the displacement of Q 
relative to Q, the medium in the subsystem S will respond to the perturbation and will relax to a new state of thermal 
equilibrium after a certain relaxation time, iroiax('S'). For example, from the results concerning the medium entropy 
and number contents in the subsystem 5* as a function of the separation between Q and Q obtained in the last few 
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sections (Fig. 1), we know that under a displacement of the relative separation of the Q and Q, the medium particles 
will exchange between S and the thermal bath S' in order to make the subsystem S under thermal equilibrium. In 
this case, the relaxation time trciax('S') corresponds to the exchange of medium particles through imaginary boundaries 
between S and the thermal bath 5". Relaxation time grows smaller as the subsystem S decreases in size [131 • For a 
small subsystem S, this relaxation time treia.x{S) can be very short. On the other hand, for heavy Q and Q in the 
subsystem, the period of Q-Q relative motion, igg, can be relatively long because of the large mass of Q and Q. The 
period tgQ can be so much greater than the medium relaxation time ireiax('S'), tgQ » trciax('S'), that the medium can 
be approximately considered as reaching a state of thermal equilibrium approximately instantaneously, at any time 
during the (supposedly slow) motion of the heavy quark Q and antiquark Q. For the medium particles, this is just 
the Born-Oppenheimer approximation for the description of the states of the medium particles in the subsystem S, as 
presented in the last few sections and used in lattice gauge calculations to obtain the medium particle configurations 
in QCD. 

In what sense can energy and entropy be considered conserved under a periodic motion of the Q and Q in the 
medium? The whole closed system consists of the subsystem S' and S, and the subsystem S is not a closed system. 
From the results of the last few sections, we know that if we move the Q closer relative to the Q in the subsystem S, then 
the medium particles (and its entropy and energy contents) will move from the subsystem S into the complementary 
part S' so as to maintain thermal equilibrium in the subsystem S. When we move the Q farther relative to the Q, 
then the medium particles, entropy, and energy contents will move from the complementary system S" back into the 
subsystem S so as to maintain thermal equilibrium in the subsystem S. For slow periodic motion of Q and Q, the 
motion can be so slow that the exchange of medium particles between S and S" can be approximated as taking place 
with no excitation of the medium. In this sense, the idealized periodic motion of Q and Q can thus be adiabatic in 
the lowest order with an "adiabatic" exchanging the energy content and the entropy content of the medium particles 
back and forth between the subsystem S and the complementary system S' . Additional interactions of the bound 
periodic Q-Q states with the medium particles that lead to the non-adiabatic excitation of both objects can then be 
considered in higher-order approximations. 

Thus, in this adiabatic picture of tgg >> t^cUxiS), the Q and Q experience the interactions from all medium 
particles which adjusts themselves (within a relaxation time which is taken to be so small as to be approximately 
instantaneous) at all instances of the dynamical motion of the Q and Q. The potential energy of the Q and Q at a 
separation r in the Hamiltonian for the Q and Q is half of the integral of the product of the local point charges of the 
Q and Q with their local potentials ^(r, R) arising from a self-consistent rearrangement of all particles when Q and 
Q are at the separation R, excluding the infinite self-energy contributions. In addition to the potential energy of the 
Q and Q, the Hamiltonian for the Q-Q system consists also of the kinetic energy of Q and Q. The Hamiltonian for 
the Q-Q system is therefore given by 



H = ^ + ^ + ^ I dr 















2mQ 


2mQ 



nr,R) 



i[e+V^(r,R)|r_R/2 + e_F(r,R)|r=R/2]. (42) 

Upon making the change of the variables from pg and pg to the center-of-mass momentum Pcm = Pq + Pg and the 
relative momentum p^j = (pg — Pg)/2, and using the solution V{r, R) of Eq. we find from the above equation 

p2 2 „2 -p.R 2 

H=—JJM +^^11 ^, (43) 

2{mQ+mQ) 2/ii.cd R 

where /Zred = '^qt^t-qU'^q + "^g) is the reduced mass. The Hamiltonian then separates into H =^ i?cM + Hn where 
Hqm is the Hamiltonian for the free motion of the composite two-body system and Hn is the Hamiltonian for the 
relative motion of Q and Q, 

From the above equation, we recognize that the potential for the Q-Q system under screening by the medium, Uqq{R), 
is given by 

^2„ ™2 

Uqq{R) = -'-^-^. (45) 



The QQ potential Uqq{R) is just the Debye screening potential plus an i?-independent constant term. 
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Based on the "adiabatic" picture of the motion of Q and Q, the Hamiltonian formulated in Eq. indeed gives 
correctly the Hamiltonian with the Debye screening potential. The i?-independent term —q^/r^ is also an important 
part of the screening contribution. We note that if we expand the Q-Q Hamiltonian for the case of small /.t representing 
the screening effects, then Eq. (|44|) becomes, 

tiR ^ = (45) 

2Mrcd R 2^i.od R 

which is the same Hamiltonian as that of the unscreened case. Thus, we reach the interesting result that in the lowest 
order of the screening parameter fi ~ l/j'j,, a properly calibrated Hamiltonian of a system under screening, with the 
shift of the level of the potential, —q^/ru, is the same Hamiltonian without screening. In practical terms, if we calculate 
the mass of a Q-Q system without screening, we expect that within the lowest order of the screening parameter /x, 
the mass eigenvalue of the system to be nearly unchanged when screening is present. Numerical calculations of the 
mass of bound L = charmonium using the potential model of Ref. [23| indeed shows that the absolute value of the 
charmonium L ~ mass changes only very slightly as a function of temperature up to 1.5Tc [soj . 

From the above discussions in the simple case of Debye screening, we observe that the potential between Q and Q, 
Uqq{R), differs from the total internal energy Utota.i{R)- Because of Eqs. (^5]) and ([57]) . they are related by 

Uqq{R) = Utotal{R) - [C/mcdium(i?) - [^medium]- (47) 

It is therefore necessary to subtract out the change of the medium internal energies [C/mcdium(-R) — C^modium] from the 
total internal energy UtotsiiR) to obtain the Q-Q potential Uqq{R) in the grand canonical ensemble. This conclusion 
for Debye screening supports a similar conclusions in the analogous lattice gauge theory, where we have proved in Eq. 
(11) of Ref. [13, 

U^qI{R,T) = U,{R,T) - [uI'Hr,T) - Ugo{T)]. (48) 

In the above equation, the superscript (1) refers to the color-singlet state of Q and Q, and [Ug^\R,T) — Uga{T)] is 
the increment of gluon energy due to the presence of Q and Q. 

From the above analysis, we conclude that the relationship of Eq. in Ref. [l^l is a rather general result for 
heavy particles under screening in the grand canonical ensemble. 

We can understand the results of Eq. (48) [or similarly (47)] from another viewpoint. In a standard description of a 
(QQ) in a medium, we simplify the dynamics by considering first the (QQ) states and the deconfined medium states 
separately as independent unperturbed states. We then include their mutual excitations as perturbative couplings. 
Thus, in the lowest-order description without perturbative couplings, the (QQ) states should be obtained without 
the excitation of the medium states of deconfined real gluons and vice versa. In the quenched approximation, the 
deconfined real gluons in the quark-gluon plasma in the Feynman diagram language are those represented by lines 

with external legs. The change of the medium internal energies, [uI^\r,T) - Ugo{T)] in Eq. (48), represents the 
excitation of the internal energy states of the deconfined real gluon medium when the separation between the Q 
and the Q changes in the grand canonical ensemble. As the unperturbed states of the Q-Q relative motion should 
be calculated without the excitation of the medium states of deconfined real gluons, we therefore need to subtract 
the change of the real gluon internal energy from the total internal energy in Eq. (48) to obtain the Q-Q potential 
U^q{R,T), when the separation between the Q and the Q changes. 

In the quenched approximation, the subtraction of this change in the internal energy of real gluons as a function 
of R does not mean that both the real and the virtual gluon degrees of freedom are frozen. Only the real gluon 
excitation energy degrees of freedom are frozen when we calculate the unperturbed Q-Q bound states of relative 
motion, for reasons we have just given. On the other hand, virtual gluons, which in the Feynman diagram language 
are represented by gluon lines with the two end points of each line joining onto other quarks and gluons, change their 
configurations as the separation R between the Q and the Q changes. These virtual gluons mediate the interaction 
between the Q and the Q. The changes in the virtual gluon configurations modify the interaction between the Q 
and the Q, resulting in the screening of the Q-Q potential. This type of virtual gluon excitation is not frozen and is 
included in the calculation of the Q-Q potential and the evaluation of bound states of the Q-Q relative motion. 

The reconfiguration of these virtual gluons can take place in either an adiabatic manner or in the opposite "diabatic" 
manner depending on whether the time scale of the relaxation of these virtual gluons is short compared to the time 
scale of the period of the Q-Q relative motion. As we explained earlier, the adiabatic description is appropriate for 
very heavy quark pairs when the relative motion of the Q and the Q is slow. This is indeed confirmed by the successful 
identification of the lattice free energy for a static Q-Q pair as the heavy quark Q-Q potential, U^q{R) at T = 
[Bali et al. Phys. Rev. D56 2566 (1997)]. For this case of T = 0, there are no free gluons nor energy excitations of 
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free gluon medium states when the separation R between the Q and the Q changes, and the total free energy of the 
system is equal to the total internal energy of the system. The real gluon energy degrees of freedom are absent and 
frozen but the virtual gluon degrees of freedom adjust themselves as R changes. 

The question whether an adiabatic or a diabatic picture is a more appropriate description for the potential arises 
also in the NN and meson-meson problems at T = in lattice gauge theory. However aside from the question of 
adiabaticity, the Q-Q, NN, and meson-meson potentials differ in their different degrees of freedom and the methods of 
calculations. The Q-Q potential studied here is a (two-body )-plus-(deconfined medium) problem, while the NN and 
meson-meson potentials at T = are (six-body )-plus- (virtual gluons) and (four-body)-plus-(virtual gluons) problems 
respectively. At T 0, a wave function treatment of the lattice gauge correlator results in the correct repulsive 
potential for the NN potential at short distances [11], while an "adiabatic" potential treatment without using the 
lattice wave function gives flat NN and meson-meson potentials [s^. On the contrary, however, another "adiabatic" 
lattice gauge meson-meson potential calculation at T = gives re pul sive and attractive inner cores when different 
internal degrees of freedom of the light quarks are taken into account [57| . Furthermore, the lattice gauge wave function 
method of [55| may need additional justifications as questions have been raised in Appendix A of 57 1 concerning its 
lattice wave function assumption. While the work of [s^ appears to give a correct description, much work remains 
to be carried out to sort out the differences of the lattice gauge calculations of [s^ , [IHl , and [131 . It remains another 
separate additional question how one can obtain definitive conclusions on the adiabaticity or diabaticity of the (two- 
body)-plus-(deconfined medium) potential at T > Tc from these (six-body)- and (four-body )-plus- (virtual gluon) 
problems at T = 0. As many unanswered questions remains to be resolved, the results of [5^ cannot yet be used, for 
the present time at least, to draw conclusions on the adiabaticity or diabaticity of the Q-Q potential examined here. 
Nevertheless, the exploration of the relationship between adiabaticity and the shape of the relative wave function is 
an interesting subject for future investigations. 



VII. AN APPROXIMATE METHOD TO SEPARATE OUT THE Q-Q POTENTIAL FROM Ui 

From the simple model of Debye screening, we observe that up to the first order of eV/T, the internal energy of 
the medium does not change, but up to the second order the internal energy increases with an increasing separation 
between Q and Q. This increase arises from the fact that the thermal equilibrium attained through the contact with 
a thermal bath in a grand canonical ensemble constrains the occupation numbers of the medium particles, and this 
newly re-arranged distribution leads to an increase in the number, the entropy, and the internal energy of the medium, 
as a function of increasing R. 

Returning now to QCD lattice gauge calculations and noting its similarities with Debye screening of Coulomb 
charges, we should therefore expect that the number, the entropy, and the internal energy of the gluon medium 
should likewise increase as function of increasing R between Q and Q. Indeed, as shown for lattice gauge calculations 
at a fixed temperature in Fig. 1(a), there is an increase in the entropy of the system as R increases, similar to the 
analogous Debye screening case shown in Fig. l{b). 

Having understood the behavior of various thermodynamic quantities, we wish to extract the Q-Q potential from 
lattice gauge results. The most reliable way is to carry out additional lattice gauge calculations to obtain Ug^\R) 
and Ugo- The Q-Q potential is then the difference of Ui{R) and Ug^\R) — Ugo, as given by Eq. (gS]). As U^g^\R) and 
Ugo in lattice gauge calculations are not yet available, we will try to use another piece of lattice gauge data to obtain 
the Q-Q potential, as least approximately. 

We note that in the Debye screening case ?7modium(-R) — t^mcdium proportional to TS{R), and in the lattice gauge 
calculations the quantity TSi{R,T) has been calculated. We can look for a similar relationship between the gluon 
internal energy and the gluon entropy for the quark-gluon plasma. If we succeed in r elating [/^^' (R, T) - Uga to 
TSi{R,T), then the Q-Q potential, can be determined from Ui{R,T) by subtraction using Eq. 

The subtraction can be carried out by noting that locally the quark-gluon plasma internal energy density e is related 
to its pressure p and entropy density a by the First Law of Thermodynamics, 

e = Ta-p, (49) 

and the quark-gluon plasma pressure p is also related to the plasma energy density e by the equation of state p{e) 
that is presumed known by another lattice gauge calculation. Thus, by expressing p as (3p/e)(e/3) with the ratio 
a(T) = 3p/e given by the known equation of state, the plasma internal energy density e is related to the entropy 
density Tcr by 

3 -I- a(T) 



14 



This is just 

where (To is the entropy density in the absence of Q and Q. Noting that the entropy of the medium for the color-singlet 
Q-Q pair is TSi ~ J drT{a — ctq) and Ugo is related to / dr Tao, the above equation leads to 

d[U^g'\R,T) - UgojT)] _ 3 TdS,{R,T) 

dV 3 + a{T) dV ' ^ ^ 

and the plasma internal energy integrated over the volume is given by 

uj^'\R,T)- Ugo {T) = 3 ^ ) TS^ {R, T) . (53) 

But TSi{R,T) has already been obtained as Ui{R,T) — Fi{R,T). The plasma internal energy is therefore equal to 

3 

'i{T) 

,(1) 
QQ 

combination of Fi and Ui given by [27[, 



f/W(i?,r) _ Ugo - T-T7^Pi{R,T) - Fi(i?,r)]. (54) 
The Q-Q potential, t/ij^A, as determined by subtracting the above plasma internal energy from C/i, is then a linear 



M^i(i?,r) ^ ul^l{R,T) - ^-A_Fi(i?,r) + 3°^J(y) ^i(^>r), (55) 

where for brevity of notation we have renamed Uqq{R,T) as Wi{R,T) and we can define the coefficient of Fi, 
fp = 3/(3 + a{T)), as the Fi fraction, and the coefficient of Ui, fu = a{T)/{3 + a{T)), as the Ui fraction. The 
potential uj^Q is approximately Fi near Tc and is approximately 3Fi/4 + J7i/4 for T > 1.5rc [l^l- 



VIII. COMPARISON OF DIFFERENT Q-Q POTENTIALS 



In the spectral function analyses, the widths of many color-singlet heavy quarkonium states broaden suddenly at 
various temperatures [2l|, [l^, ■ In the most precise calculations for J /if; using up to 128 time- like lattice slices, the 
spectrum has a sharp peak for O.TSTc I^T < 1.62Tc and a broad structure with no sharp peak for l.TOTc <T< 2.33rc 
[2ll |. The spectral peak at the bound state has the same structure and shape at 0.78Tc as it is at 1.62Tc. If one can 
infer that J/tp is stable and bound at 0.78Tc, then it would be reasonable to infer that J/ip is also bound and stable at 
1.62 Tc- The spectral function at l.TOTc has the same structure and shape as the spectral function at 2.33rc- If one can 
infer that J/ip is unbound at 2.33rc, then it would be reasonable to infer that J/ip become already unbound at 1.70 Tc- 
We can define the spontaneous dissociation temperature of a quarkonium as the temperature at which the quarkonium 
changes from bound to unbound and dissociates spontaneously. Thus, from the shape of the spectral functions, the 
temperature at which the width of a J/iJj quarkonium broadens suddenly from 1.62Tc to l.TOTc corresponds to the 
J/zispontaneous dissociation temperature. Spontaneous dissociation temperatures for Xc and Xb have been obtained 
in [22I. [ssj . We list the heavy quarkonium spontaneous dissociation temperatures obtained from spectral analyses in 

quenched QCD in Table I. They can be used to test the potential models of Wi(= U^q), Fi, and Ui. 



Table I. Spontaneous dissociation temperatures obtained from different analyses. 





Quenched QCD 


FuU QCD (2 flavors) 


States 


Spectral Analyses 


Wi 


Fi 


Ui 


Wi 


Fi 


Ui 




1.62-1.70TJ 


1.62 


1.40 Tc 


2.60 Tc 


1.42 Tc 


1.21 Tc 


2.22 Tc 


Xc 


below l.lTc'i 


unbound 


unbound 


1.18 Tc 


1.05 Tc 


unbound 


1.17 Tc 






unbound 


unbound 


1.23 Tc 


unbound 


unbound 


1.11 Tc 






4.1 Tc 


3.5 Tc 


5.0 Tc 


3.40 Tc 


2.90 Tc 


4.18 Tc 


Xb 


1.15-1.54Tc» 


1.18 Tc 


1.10 Tc 


l.T3Tc 


1.22 Tc 


1.07 Tc 


1.61 Tc 






1.38 Tc 


1.19 Tc 


2.28 Tc 


1.18 Tc 


1.06 Tc 


1.47 Tc 



fRef.[2lj, «Ref.[22l, «Ref.[35] 
To evaluate the Q-Q potentials in quenched QCD, we use the free energy Ti and the internal energy Ui obtained by 
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Kaczmarek et al. |33j where Fi and Ui can be parametrized in terms of a screened Coulomb potential with parameters 
shown in Figs. 2 and 3 of Ref. [23|- For the ratio a{T) from the plasma equation of state in Eq. ([55]) . we use the 
quenched equation of state of Boyd et al. [sl] for quenched QCD. The quantity a{T) = Sp/e and the Ui and Fi 
fractions as a function of T are shown Fig. l{b) and Fig. 1(c) of Ref. [13] respectively. These quantities allows the 
specification of the Wi = U^q potential as a function of temperature. 

Using quark masses rric = 1.41 GeV and mt — 4.3 GeV, we can calculate the binding energies of heavy quarkonia and 
their spontaneous dissociation temperatures using different potentials in quenched QCD. As a function of temperature, 
the binding energies and wave functions of charmonia have been presented in Fig. 6 and 7 of Ref. [27j respectively. 
The bounding energies of bottomia have been presented in Figs. 8 and 9, and the wave functions in Fig. 10 of Ref. 
[131. We show the root-mean-square Q-Q separation yJJW^ of L = charmonium calculated with mc = 1.41 GeV in 
the Wi{R) potential in quenched QCD as the sohd curve in Fig. 4. As one observes, the root-mean-squared separation 
-Rrms is about 1 fm for T ~ Tc and it increases to about 4 fm at T = 1.6Tc before it becomes unbound. The large 
separation is expected for systems with a weak binding, in analogous to the halo nuclei with neutrons in weak binding 
observed in nuclear physics [s^ . There is the question whether charmonium with such a large separation between Q 
and Q may survive in QGP. The dissociation cross sections for these quarkonia by collision with gluons have been 
calculated and found to be a function of the gluon collision energy, as shown in Fig. 13 of Ref. [27j . 




1 1.2 1.4 1.6 1.8 



FIG. 4: The root-mean-squared Q-Q separation of L = charmonium as a function of T/Tc in quenched QCD (solid curve), 
and in full QCD with two favors (dashed curve). 

From the binding energy of a quarkonium as a function of temperature, one can obtain the temperature at which 
the quarkonium binding energy vanishes. This is the temperature for the spontaneous dissociation of the quarkonium, 
as the quarkonium at this temperature will dissociate spontaneously. We list in Table I the heavy quarkonium 
spontaneous dissociation temperatures calculated with the Wi potential, the Fi potential, and the Ui potential, in 
quenched QCD. 

The J/il) and Xb spontaneous dissociation temperatures obtained with the Wi potential in quenched QCD are 
found to be \.(!)2Tc and l.lSTc respectively. Spectral analyses in quenched QCD give the spontaneous dissociation 
temperature of 1.62-1.70Tc for J/ijj [2l| and 1.15-1.54Tc for Xb [la]- Thus, spontaneous dissociation temperatures 
obtained with the Wi potential agree with those from spectral function analyses. This indicates that the Wi potential, 
defined as the linear combination of Ui and Fi in Eq. (j55p . may be the appropriate Q-Q potential for studying the 
stability of heavy quarkonia in quark-gluon plasma. 



IX. Q-Q POTENTIAL FOR FULL QCD WITH TWO FLAVORS 



The interaction energy between a heavy quark and a heavy antiquark in the color-singlet state in two-flavor full 
QCD was studied by Kaczmarek and Zantow [3J|. In full QCD with 2 flavors, Fi and Ui can be represented by a 
color-Coulomb interaction at short distances and a completely screened, constant, potential at large distances as given 
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in Rcf . [28| , although other alternative representations have also been presented [36|, |37[ . The transitional behavior 
linking the two different spatial regions can be described by a radius parameter ro{T) and a diffuseness parameter 
d{T), as in the Wood-Saxon shape potential in nuclear physics, 



4 (T (T) 

{F,, U,}{R, T) = --^y /(i?, T) + C{T)[l - /(i?, T)], 



3 R 



where 



exp{(i?-ro(T))MT)} + l- 



(56) 



(57) 



In principle, it is necessary to specify only the temperature dependence of Fi{R, T) as the internal e nerg y Ui{R, T) can 
be obtained from Fi and its derivative with respect to T. In practice, as Kaczmarek and Zantow [34| have obtained 
Ui (i?, T) by a careful numerical differentiation, it is convenient to parametrize the internal energy in the above simple 
form for practical calculations. 




2.0 2.5 3.0 
T/T 



FIG. 5: The parameters C, ro, and d for the color-singlet free energy Fi{R,T) in two-flavor QCD as given in Eq. (|56|l . 



In searching for the coupling constant that fits the lattice quantities, we found that the value of as centers 
around 0.3. The fit to the lattice gauge quantities does not change significantly whether we allow as to vary. It is 
convenient to keep the value of as to be 0.3 so that there are only three parameters for each temperature. 

For the free energy Fi{R,T) in two-flavor QCD [33], the set of parameters C, ro, and d are shown in Fig. 5, and 
the corresponding fits to Fi are shown in Fig. 6. For the internal energy Ui{R,T) in two-flavor QCD the set of 
parameters C, ro, and d are shown in Fig. 7, and the corresponding fits to the lattice gauge internal energy Ui results 
shown in Fig. 8. 

If the thermodynamical quantity Fi or Ui are treated as a potential, then the quantity C{T) is an approximate 
measure of the depth of the potential measured from the flat potential surface at large distances relative to the 
potential well at short distances. For the free energy Fi, the C{T) parameter has the value of about 1 GeV at 
T ~ O.STc, and it decreases to 0.5 GeV at Tc. The free energy as a potential will has a well depth of about 0.5 GeV 
for T close to Tc and the well depth decreases to about 0.1 GeV at T 2Tc. 

One notes that there is a significant change in the slopes of C(T) at T ~ Tc for Fi. As a consequence, the parameter 
C(T) for Ui exhibits a peaks at T ~ Tc. The transitional radius ro for Fi decreases gradually from about 0.6 fm 
to about 0.15 fm and the diffuseness parameter d decreases slowly from 0.3 fm to about 0.15 fm, as temperatures 
decreases from 0.7Tc to 4Tc. 
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T/T 




R (fm) 



FIG. 6: The symbols represent the color-singlet free energy, Fi{R, T), for two-flavor QCD [sj, and the curves are the fits using 
the screened potential, Eq. (|56|) . with parameters given in Fig. 5. 




FIG. 7: The parameters C, ro, and d for the color-singlet internal energy (7i(R, T) as given in Eq. (|56p . 



For the internal energy Ui, the parameter C{T) is quite large, attaining the value of about 3 GeV for T close 
to Tc. This indicates that if Ui is used as a potential, the potential depth at temperatures close to the transition 
temperature is of order 3 GeV, which is a very deep potential indeed. The parameter C(T) decreases to about 0.8 
GcV when the temperature exceeds about 1.5 Tc. The transition radius vq is about 1 fm for T close to 0.8 Tc, and it 
decreases to about 0.2 fm for T ~ 4Tc. The diffuseness parameter d{T) for the internal energy decreases substantially 
at temperatures below Tc, but maintains a relatively constant values of 0.1 to 0.2 fm for T greater than Tc- 

The comparison in Fig. 6 and 8 shows that the free energy Fi and the internal energy Ui with the set of param- 
eters in Figs. 4 and 6, adequately describe the lattice-gauge data and can be used to calculate the eigenvalues and 
eigenfunctions of heavy quarkonia. 

In our description, the Q-Q potential is a linear combination of Ui and Fi with coefficients dependin g o n the 
equation of state. The equation of state in full QCD with two flavors has been obtained by Karsch et al. [47|. We 
show their results of e/T^ and Sp/T"^ in Fig. 9(a). The ratio a(T) = 3p/e and the Ui and Fi fractions as a function 
of T are shown in Figs. 8{b) and 8(c) respectively. Similar to the case of quenched QCD, the Fi fraction is close to 
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unity near Tc and it decreases to 3/4 at large temperatures while Ui fraction is nearly zero at Tc and it increases to 
about 1/4 at high temperatures. These quantities, together with Fi and Ui specifies the Q-Q potential for bound 
state calculations. 



X. HEAVY QUARKONIA IN QUARK-GLUON PLASMA 

In full QCD with two flavors, the transition temperature is Tc = 202 MeV [U. To calculate the charmonium energy 
levels, we employ a quark mass mc = 1.41 and rUb = 4.3 GeV. 

Energy levels of charmonium states calculated with different potentials in full QCD with two flavors arc shown in 
Fig. 10 as a function of the temperature in units of Tc. The J and rjc states are weakly bound and they dissociate at 
1.21 Tc in the Fi potential, at 1.42Tc in the Wi potential, and at 2.22Tc in the Ui potential. The Xc state dissociates 
below Tc in the Fi potential, at 1.05Tc in the Wi potential, and at l.YITc in the Ui potential. At temperatures 
slightly greater than Tc, they are weakly bound in the Fi potential but are strongly bound in the Ui potential, with 
a binding energy of about 0.7 GeV at l.lTc. The binding of the states in the Wi potential lies between these two 
limits. 

The dissociation temperature for J and rjc in full QCD with two flavors have been examined in the spectral 
function analysis by Aarts et al. j60j . As the calculations have been carried out only with a small lattice volume 
and a small set of statistics, there are possible systematic uncertainties which prevented a precise determination of 
the pseudocritical temperature. Preliminary results indicate that the J/ij} state may be bound up to about 2.Tc [60| . 
More definitive results will await a greater lattice volume and larger statistics. 

Wc also carry out the analysis of bottomium Is, 2p, and 2s states. Fig. 11 gives the state energies as a function of 
T/Tc for different potentials. The eigenenergies of T and rji, in two-flavor QCD are about -0.1 GeV at T = 1.1 Tc in 
the Fi potential, and are about -1.0 GeV in the Ui potential. The eigenenergies in the Wi lie in between those of the 
Fi and Ui potentials. These states dissociate spontaneously at 2.9rc in the Fi potential, 3.40 in the Wi potential, 
and about 4 to 5Tc in the Ui potential. 

In full QCD with two flavors, the Xb state dissociates at l.OTTc in the Fi potential, at 1.22 in the Wi potential, and 
at 1.61 in the Ui potential. The T' and 77^ state dissociates at 1.06Tc in the Fi potential, at 1.18 in the Wi potential, 
and at 1.47 in the Ui potential. 
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FIG. 9: The equation of state in full QCD with two flavors: (a) the quantity e/T* and 3p/T*, (b) the ratio a(T) = 3p/e, and 
(c) the U\ and F\ fractions in the Q-Q potential C/^'i- 
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FIG. 10: Energy levels of charmonium in the quark-gluon plasma as a function of temperature calculated with the Fi{R,T), 
Wi{R,T), and Ui{R,T) potentials in two-flavor QCD. Fig. 10(a) is for the J/tp and rjc state, Fig. 10(&) is for the Xc state, and 
Fig. 10(c) is for the ip' state. 



In Table I we list the dissociation temperatures of different quarkonia obtained in full QCD with two flavors. A 
comparison of the dissociation temperatures from the quenched QCD and full QCD with two flavors show the effects 
of the dynamical quarks. Dynamical quarks increases the degree of screening, but at the same time, lowers the phase 
transition temperature. They lead to a more diffused potential with a greater screening length. As a consequence, 
the binding energy of the Is states is lowered and the dissociation temperature decreases. We can chooses the Wi 
potential as the more appropriate potential for the Q-Q pair, as it gives the best agreement with spectral function 
analysis in quenched QCD. For this Wi potential, the dissociation temperature decreases from 1.62Tc in quenched 
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FIG. 11: Energy levels ol bottomium in the quark-gluon plasma as a function of the temperature calculated with the -Fi(R, T), 
VFi(R, T), and Ui{R,T) potentials in full QCD with two flavors. Fig. 11(a) is for the T and rjt state, Fig. 11(6) is for the xt 
state, and Fig. 11(c) is for the T' state. 



QCD to lA2Tc in full QCD with two flavors. The effects of the dynamical quark in full QCD leads to a slightly weaker 
binding for J/ip in the plasma. For the x states, the effects of the the additional quark screening does not modify the 
dissociation temperature substantially. The additional screening tends to move the centrifugal barrier for the / = 1 
state to a smaller radial distance with a slightly higher barrier, resulting in a very slight increase in the dissociation 
temperature. 



XI. QUARK-DRIP LINES IN QUARK-GLUON PLASMA 

To examine the stability of a color-singlet Q-Q pair, we consider the quark mass mg as a variable and evaluate 
the spontaneous dissociation temperature as a function of the reduced mass /Xred = t^Q'niQ I {"niq + ?tiq). The quark 
drip lines calculated with the Fi, VFi, and Ui potentials in quenched QCD are shown in Fig. 12. A state is bound 
in the (T /Tc, /.iicd) space above a drip line and is unbound below the drip line. Spectral function analysis gives the 
spontaneous dissociation temperature of 1.62-1.70Tc for J/-0 [2ll. [2^ and 1.15-1.54Tc for Xb [13, |3^. If one takes the 
charm quark mass to be 1.41 GeV and the bottom quark mass to be 4.3 GeV, the spectral function results can be 
represented by the solid-circle symbols in Fig. 12. They fall on the drip line curves obtained with the Wi potential, 
indicating that the Wi is the appropriate Q-Q potential to use for bound state problems. 

In quenched QCD, however, the quark-gluon plasma is assumed to consist of gluons only. As dynamical quarks 
may provide additional screening, it is necessary to consider the case with dynamical quarks. Accordingly, we use the 
Fi, Wi, and Ui potentials evaluated in full QCD with 2 flavors [13, E^j to determine the drip lines in Fig. 13. The 
drip lines for the Ui potential lies lower than that of the Wi , which in turn lies lower than the drip line of the Fi 
potential. In comparison with quenched QCD results, the Is drip line in full QCD is shifted to lower temperatures 
while the Ip drip line in full QCD is only slightly modified. 

We shall use the results from the Wi potential in full QCD with two flavors to discuss the question of quarkonium 
stability in quark-gluon plasma, as the Wi potential has been found to give results in agreement with spectral function 
analyses in quenched QCD. For heavy quarkonia, results in Table I and Fig. 13 obtained with the Wi potential in 
fifll QCD with two flavors indicate that J/tp, Xc, , Xb, and T' may be boimd in the plasma up to 1.42Tc, l.OSTc, 
3.40Tc, 1.22rc, and l.lSTc respectively. 

The variation of the drip lines with the reduced mass allows us to examine the stability of quarkonia containing 
quarks of various masses. We need to know the effective masses of different quarks in the quark-gluon plasma. Due 
to its strong interaction with other constituents, a light quark becomes a dressed quasiparticle and acquires a large 
quasiparticle mass. In the low temperature region where the spontaneous chiral symmetry breaking occurs with 
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FIG. 12: The quark drip lines in quenched QCD calculated with the Fi, Wi, and Ui potentials. The symbols represent results 
from lattice gauge spectral function analyses. 




FIG. 13: The quark drip lines in 2-flavor QCD calculated with the Wi, Fi and Ui potentials. 



{tptfj) ^ 0, the quasiparticle mass is ^ (current quark mass)], where g is the strong couphng constant 

and (i/'V') the quark condensate [H, 0, [621. This quasiparticle mass is the origin of the constituent-quark mass in 
non-relativistic constituent quark models[43l. [H, [6l, [gj]- In the high temperature perturbative QCD region, the 
quasiparticle mass is rriq ^ gT/^/6, which is of the order of a few hundred MeV [65l |. 

As the restoration of chiral symmetry is a second order transition, (V"^) decreases gradually as the temperature 
increases beyond Tc. The light quark quasiparticle mass associated with (^/'f/') will likewise decrease gradually from 
the constituent-quark mass value to the current-quark mass value when the temperature increases beyond Tc. This 
tendency for the quasiparticle mass to decrease will be counterbalanced by the opposite tendency for the quasiparticle 
'thermal mass' to increase with increasing temperature. As a result of these two counterbalancing tendencies in the 
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region of our interest, < T < 2Ti,, the effective mass of the hght quarks are relatively constant. By examining the 
effects of the light quark quasiparticle masses on the quark-gluon plasma equation of state, Lcvai et al. (66| . Szabo 
et al. [13, and Ivanov et al. [bI] estimate that niq is about 0.3 to 0.4 GeV at Tc < T < 2Tc. As in the case of T = 0, 
where light quarks with a constituent-quark mass of about 350 MeV mimic the effects of chiral symmetry breaking 
and non-relativistic constituent quark models have been successfully used for light hadron spectroscopy [H, [H] , so 
the large value of the estimated quasiparticle mass (from 0.3 to 0.4 GeV) may allow the use of a non-relativistic 
potential model as an effective tool to estimate the stability of light quarkonia at Tc < T < 2Tc. It will be of interest 
to investigate the relativistic effects [4^, |4^ in the future. 

For light quark masses of 0.3 to 0.4 GeV, we can estimate from the results in Fig. 13 for the Wi potential that 
as a quarkonium with light quarks has a reduced mass of 0.15-0.2 GeV, it may be bound at temperatures below 
(1.05 — 1.07)Tc. An open heavy quarkonium with a light quark and a heavy antiquark or a light antiquark and a 
heavy quark have a reduced mass of about 0.3-0.4 GeV and may be bound at temperatures below (1.11 — 1.19)Tc. 

Another lattice gauge calculation gives mq/T = 3.9 ± 0.2 at 1.5Tc [i^, which implies that at T=:1.5Tc (or about 
0.3 GeV), the quark mass will be ~1.2 GeV for {u,d,s) quarks. Such a 'hght' quark quasiparticle mass appears to 
be quite large and may be uncertain, as the plasma will have approximately equal abundances of 'light' and charm 
quarks, which is however not observed. There may also be difficulties in reproducing the plasma equation of state. 
With this mass, a 'light' quarkonium will have a reduced mass of 0.6 GeV, and the quarkonium may be bound at 
temperatures below '--^1.3irc- 

In either case, the drip lines of Fig. 13 for full QCD with 2 flavors obtained with the Wi potential do not support 
bound QQ states with light quarks beyond 1.5Tc. A recent study of baryon-strangeness correlations suggests that the 
quark-gluon plasma contains essentially no bound QQ component at 1.5Tc [i^ . 

XII. CONCLUSIONS AND DISCUSSIONS 

The degree to which the constituents of a quark-gluon plasma (QGP) can combine to form composite entities is an 
important property of the plasma. To study the composite nature of the plasma, we need to examine the stability 
of quarkonium in quark-gluon plasma which depends on the Q-Q potential. We seek to extract the Q-Q potential 
from thermodynamical quantities obtained in lattice gauge calculations. For such a purpose, we need the relationship 
between the Q-Q potential and the internal energy obtained in lattice gauge calculations. Such a relationship was 
derived previously in Ref . p7l | . We would like to gain additional support by examining whether a similar relationship 
exists between the Q-Q potential and the internal energy in an analogous, but not identical, case of Debye screening. 

We find that in adiabatic motion of Q and Q under Debye screening, (1) the potential for the Q and Q in the 
Schrodinger equation contains the interactions that act on Q and Q, (2) this Q-Q potential under Debye screening is 
only part of the total internal energy of the system, (3) the other part of the internal energy is the internal energy of 
the medium particles, and (4) many thermodynamical quantities such as the number, the entropy, and the internal 
energy of the medium particles increases with the separation between Q and Q in the grand canonical ensemble. 
Therefore, to obtain the Debye screening potential between two static charges, it is necessary to subtract out the 
internal energy of the medium particles from the total internal energy. These results supports a similar conclusion 
reached earlier in the analogous lattice gauge theory [23| . 

We are thus led to obtain the Q-Q potential in the quark-gluon plasma by subtracting out the internal energy 
of the medium particles from the total internal energy in the grand canonical ensemble. We proposed a method 
to subtracting out the internal energy of the medium by making use of the equation of state of the quark-gluon 
plasma obtained in an independent lattice gauge calculation [23|. The potential can then be represented as a linear 
combination of Ui and Fi, with coefficients depending on the quark-gluon plasma equation of state. The proposed 
potential in the quenched approximation is found to give dissociation temperatures that agree with those from spectral 
function analyses. It can be generalized to the case of full QCD to discuss quarkonium states in the plasma. 

The knowledge of the single-particle states using potentials extracted from lattice gauge calculations in full QCD 
can then be used to examine the limit of stability of both heavy and light quarkonia and to determine the location of 
the quark drip lines. 

The quark drip lines allows one to ascertain the degree of stability of heavy and light quarkonia when the masses 
of the quarks are known. J ftp, Xc Xb, and T' are found to be stable in the plasma and dissociate at different 
temperatures. The characteristics of the quark drip lines severely limit the region of possible quarkonium states with 
light quarks to temperatures close to the phase transition temperature. Various estimates give a light quark mass of 
about 0.3-0.4 GeV |66l. IbtI. [68j . which is not very different from the constituent quark masses in non-relativistic quark 
models of hadrons. Bound quarkonia with light quarks may exist very near the phase transition temperature if their 
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effective quark mass is of the order of 300-400 McV and higher. 

The author thanks Drs. H. Crater and Su-Houng Lee for helpful discussions. This research was supported in part 
by the Division of Nuclear Physics, U.S. Department of Energy, under Contract No. DE-AC05-00OR22725, managed 
by UT-Battle, LLC and by the National Science Foundation under contract NSF-Phy-0244786 at the University of 
Tennessee. 



K. Adcox et al, PHENIX Collaboration, Formation of Dense Partonic Matter in Relativistic Nucleus-Nucleus Collisions 
at RHIC: Experimental Evaluation by the PHENIX Collaboration, nucl-cx/0410003 

J. Adams et al., STAR Collaboration, Experimental and Theoretical Challenges in the Search for the Quark Gluon Plasma: 
The STAR Collaboration's Critical Assessment of the Evidence from RHIC Collisions, |nucl-ex/0 501009 
L Arsene et al., BRAHMS Collaboration, Quark Gluon Plasma an Color Glass Condensate at RHIC? The perspective from 
the BRAHMS experiment, |nucl-ex/0410020 , 

B. B. Back et al., Phobos Collaboration, The PHOBOS Perspective on Discoveries at 7?ff/C, nu cl-ex/0410022[ 
K. P. Das and R. C. Hwa, Phys. Lett. 68B, 459 (1977); R. C. Hwa, Phys. Rev. D 22, 1593 (1980). 

R. C. Hwa and C.B. Yang, Phys. Rev. C70, 024904 (2004). 

V. Greco, C. M. Ko, and P. Levai, Phys. Rev. Lett. 90, 202302 (2003). 

V. Greco et al., Phys. Rev. C 68, 034904 (2003). 

R. J. Fries, B. MuUer, C. Nonaka, and S. A. Bass, Phys. Rev. Lett. 90, 202303 (2003); Phys. Rev. C 68, 044902 (2003). 

R. J. Fries, S. A. Bass, and B. Miiller, Phys. Rev. Lett. 94, 122301 (2005). 

P. Braun-Munzinger, J. Stachel, J. P. Wessels, and N. Xu, Phys. Lett. B 344, 43 (1995). 

P. Braun-Munzinger, J. Stachel, J. P. Wessels, and N. Xu, Phys. Lett. B365, 1 (1996). 

J. Cleymans and K. Redlich, Nucl.Phys. A661, 379 (1999). 

P. Braun-Munzinger and J. Stachel, Nucl. Phys. A606, 320 (1996). 

J. Cleymans and K. Redlich, Phys. Rev. C60, 054908 (1999). 

J. Cleymans, H. Oeschler, and K. Redlich, J. Phys. G25, 281 (1999). 

P. Braun-Munzinger and J. Stachel, J. Phys. G28, 1971 (2002). 

P. Braun-Munzinger, J. Stachel, and C. Wetterich, Phys. Lett. B596, 61 (2004). 

P. Braun-Munzinger, K. Redlich, and J. Stachel, flucl-th / 03040 1 3 , in Quark Gluon Plasma 3, (eds.) R. Hwa and X.N. 
Wang. 

J. Cleymans, H. Oeschler, K. Redlich, and S. Wheaton, Phys. Lett. B615 , 50 (2005). 

M. Asakawa, T. Hatsuda, and Y. Nakahara Nucl. Phys. A715, 863 (2003); M. Asakawa and T. Hatsuda, Phys. Rev. Lett. 
92, 012001 (2004); T. Hatsuda, hep-ph/0509306 

S. Datta, F. Karsch, P. Petreczky and I. Wetzorke, Phys. Rev. D69, 094507 (2004) and J. Phys. G31, S351 (2005). 

E. V. Shuryak and L Zahed, Phys. Rev. C70, 021901 (R) (2004); Phys. Rev. D70, 054507 (2004); E. V. Shuryak, Nucl. 

Phys. A750, 64 (2005). 

G. E. Brown, C. H. Lee, M. Rho, and E. Shuryak, Nucl. Phys. A 740, 171 (2004). 
T. Matsui and H. Satz, Phys. Lett. B178, 416 (1986). 

C. Y. Wong, Phys. Rev. C 65, 034902 (2002); J. Phys. G28, 2349 (2002). 
C. Y. Wong, Phys. Rev. C72, 034906 (2005). 

C. Y. Wong, in Proceedings of XVHth Recontre de Blois, Chateau de Blois, France, May 15-20, 2005, edited by M. 



Haguenauer, B. Nicolescu, and J. T. Thanh Van, p. 467, hep-ph/0509088 
C. Y. Wong, J. Phys. G32, S301 (2006). 

C. Y. Wong and H. W. Crater, Phys. Rev. D75, 034505 (2007). 

S. Digal, P. Petreczky, and H. Satz, Phys. Lett. B514, 57 (2001); Phys. Rev. D64, 094015 (2001). 

O. Kaczmarek, F. Karsch, and P. Petreczky, and F. Zantow, Phys. Lett. B543, 41 (2002); F. Zantow, O. Kaczmarek, F. 
Karsch, and P. Petreczky, hep-lat/0301015 

O. Kaczmarek, F. Karsch, P. Petreczky, and F. Zantow hep-lat/030912l| 



O. Kaczmarek and F. Zantow, Phys. Rev. D71, 114510 (2005), and O. Kaczmarek and F. Zantow, [liep-lat/0506019 



K. Petrov, Eur. Phys. J. C43, 67 (2005); A. Mocsy, Talk presented at Quark Matter Conference, 2005, |hep-ph/0510135l 
S. Digal, O. Kaczmarek, F. Karsch, and H. Satz, Eur. Phys. J. C43, 71 (2005). 
W.M. Alberico, A. Beraudo, A. De Pace, and A. Molinari, hep-ph/0507084 
D. Blaschke, O. Kaczmarek, E. Laermann, and V. Yudichev, Eur. Phys. J. C43, 81 (2005). 
T. Song and S. H. Lee, Phys. Rev. D 72, 034002 (2005). 

C. DeTar, Phys. Rev. D 32, 276 (1985); C. DeTar, Phys. Rev. D 37, 2328 (1988). 
T. H. Hansson, S. H. Lee, and L Zahed Phys. Rev. D 37, 2672 (1988). 

Yu. A. Simonov, Phys. Atom. Nucl. 58, 309 (1995). Yad. Fiz. 58N2, 357 (1995); Yu.A.Simonov, Phys. Lett. B619 293 
(2005). 

T. Hatsuda and T. Kunihiro, Phys. Rev. Lett. 55, 158 (1985). 

V. Koch, A. Majumder, and J. Randrup, Phys. Rev. Lett. 95, 182301 (2005). 

F. Karsch, S. Ejiri, and K. Redlich, Nucl. Phys. A774, 619 (2006). 



24 



F. G. Werner and J. A. Wheeler, Phys. Rev. 106, 126 (1958). 

F. Karsch, E. Laermann, and A. Peikert, Phys. Lett. B478, 447 (2000); F. Karsch, Nucl.Phys. A698, 199 (2002). 
H. W. Crater and P. Van Alstine, Phys. Rev. D70, 034026, (2004). 

H. W. Crater, C. Y. Wong, and P. Van Alstine, Phys. Rev. D74, 054028 (2006). 
A. Mocsy and P. Petreczky, Phys. Rev. D73, 074007 (2006). 
A. Mocsy and P. Petreczky, J. Phys. G32, S515 (2006). 

L. D. Landau, Izv. Akad. Nauk. SSSR. ser. fiz. 17, 51 (1953); L. D. Landau and S. Z. Belenkii, Usp. Fiz. Nauk 56, 309 
(1955). 

L. D. Landau and L. M. Lifshitz, Statistical Physics, Pergamon Press, Oxford, 1959, Chapter I. 

G. Bah, K. Schilling, and A. Wachter, Phys. Rev. D56 2566 (1997). 
N. Ishii, S. Aoki, and T. Hatsuda, hep-lat/0610002 



T. Takahashi, T. Doi, and H. Suganuma, hep-lat/0601006 and T. Doi, T. Takahashi, and H. Suganuma, |hep-lat /0601008 
W. Detmold, K. Orginos, and M. J. Savage, jhep-lat /0703009| 
G. Boyd et al., Nucl. Phys. B 469, 419 (1996) 



L Tanihata, Prog. Part. Nucl. Phys. 35, 505 (1995); V. Rotival and T. Duguet, 'nucl-th/0702050' 

G. Aarts, C. AUton, R. Morrin, A. O Cais, M. B. Oktay, M. Peardon, J.-L SkuUerud, Nucl.Phys.Proc.Suppl. 153 296, 



(2006) [ hep-lat/0511028 



Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961), and 124, 246 (1961). 
A. P. Szczepaniakl and E. S. Swanson, Phys. Rev. Lett. 87, 072001 (2001). 
T. Barnes and E. S. Swanson, D46, 131 (1992). 

C. Y. Wong, T. Barnes and E. S. Swanson, Phys. Rev. C65, 014903 (2002). 

H. A. Weldon, Phys. Rev. D26, 2789 (1982); K. Kajantie, M. Laine, K. Rummukainen, and Y. Schroder, Phys. Rev. D67, 
105008 (2003); J. P. Blaizot, E. lancu, and A. Rebhan, Phys. Rev. D68, 025011 (2003). 
P. Levai and U. Heinz, Phys. Rev. C57, 1879 (1998); 
K. K. Szabo and L Toth, JHEP 0306 008 (2003). 

Yu. B. Ivanov, V.V. Skokov, and Toneev, Phys. Rev. D71, 014005 (2005). 

P. Petreczky, F. Karsch, E.Laemann, S. Stickan, and L Wetzorke, Nucl. Phys. B, Proc. Suppl. 106, 513 (2002). 



